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Abstract 

It is well known that the squeezing spectrum of the field exiting a nonlinear cavity can be directly 
obtained from the fiuctuation spectrum of normally ordered products of creation and annihilation 
operators of the cavity mode. In this article we show that the output field squeezing spectrum 
can be derived also by combining the fiuctuation spectra of any pair of s~ordered products of 
creation and annihilation operators. The interesting result is that the spectrum obtained in this 
way from the linearized Langevin equations is exact, and this occurs in spite of the fact that no 
s-ordered quasiprobability distribution verifies a true Fokker-Planck equation, i.e., the Langevin 
equations used for deriving the squeezing spectrum are not exact. The (linearized) intracavity 
squeezing obtained from any s-ordered distribution is also exact. These results are exemplified in 
the problem of dispersive optical bistability. 
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I. INTRODUCTION 



An appropriate tool for studying the fluctuations of a quantum light field E (t) 
is their spectrum, which is defined as the Fourier transform of the field correlations 
(^E (t + t) ,E (t)^, where {U, V) = {UV) — {U) {V). Here we are concerned with the study 
of the amount of squeezing provided by nonlinear cavities, i.e., optical cavities containing a 
nonlinear medium and pumped by some input field. Nonlinear cavities are known to produce 
large amounts of squeezed light for a particular frequency or band of frequencies in the field 
exiting the cavity, as a result of the interference at the cavify output mirror between the 
pa^iall, .<,ue.ed .ode and t.e ..e™. ,.od. 

When calculating the fluctuations of the output field quadratures, one uses the squeezing 
spectrum 0, Q] defined as 



-oo 



I dre-^-^ (: X^* (t + r) , (t) :) , (1) 

— oo 

where the field quadrature 

^out _ Q-out jt) e ^'^ + Q-out jt) e"^ ^2"^ 

being dl^^ and dout the creation and annihilation operators of the output field that verify 
dout {t) , diut {t') = ^ {t ~ t'), and if an arbitrary phase. In Eq. ((H) the label "out" refers 
to the field exiting the nonlinear cavity, : : denotes normal and time ordering, and the term 
i corresponds to the shot noise level. 



A. Squeezing spectrum from the generalized P distribution 

The use of Eq. requires relating the correlations of the field outside the cavity -which 
are the ones that are actually detected- with those of the intracavity field, which are readily 
calculated by solving the master equation for the particular system. This relation is given 
by the input-output theory (5|], which, when the nonlinear cavity is fed with a coherent or 
vacuum field, states that 6| 

(: flout (t + r), dout (t) :) = 7^^^ (: a (t + r) , a (t) :) , (3) 
(: flLt {t + ^), flout (t) :^ = 7„^t (: (t + r) , a (t) :) , (4) 
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where 7o^t represents the cavity loss rate of the field intensity at the output mirror, and 
and d are the creation and annihilation operators of the intracavity field, which verify 
[d(t) , d^ (t)] = 1. We recall that no simple relations like Eq. (j3l4p exist that relate corre- 
lations between outgoing fields and intracavity fields unless those are calculated in normal 
order. 

The dynamics of the nonlinear cavity is usually described, although not necessarily, 
through Langevin equations obtained from a Fokker-Planck equation. Then the form of 
Eqs. fj3l4p suggests the use of the generalized P representation 8|, which we denote by P, 
as in this case 

(: a (t + r) , d (t) :) = {a (t + t) , a (t))^ , (5) 
(: dt (t + r) , d (t) :) = (/? (t + r) , a (t))^ , (6) 

where a and (3 are independent c-numbers associated to d and d^ respectively that verify 
= {a)p, and {f)p denotes the average value of any function / {a, (3) calculated in the 
V representation. Then, making use of the above equations one can write the squeezing 
spectrum as 

5rM = ^ + 7outVp(^,^) , (7) 

+ 00 

— oo 

with X^ (t) = [a (t) e~"^ + f3 (t) e"^] /2, which is a well known and widely used result. 

B. Squeezing spectrum from other quasiprobability distributions 

The above presentation suggests that the use of the V distribution is mandatory in 
order to derive the squeezing properties of the output field. (The usual Glauber-Sudarshan 
P distribution presents well known problems that will be recalled below.) Nevertheless, 

Q 

and this is at the heart of our work, a suitable combination of two s-ordered [9|] two-time 
correlations, for example antinormally (s = —1) and symmetrically (s = 0) ordered as 
obtained by using the Q (Husimi) and W (Wigner) representations 3|, |5|,|l0|, leads to the 
same result. 
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Consider the obvious property 



2 (a (t + r) , a (t))^ = (a (t + r) , a (t))^ + (a (t + r) , a (t))^ 



(8) 



where the subscript indicates the quasiprobabihty distribution used for 
obtaining the correlations. Then it follows that (: d (t + r) , a (t) :) = 

2{a{t + T),ait))^ - {a{t + T),ait))Q, and then (: X^^' (t + r) , X^^^t) :] 



7, 



out 



2 {X^ (t + r), X^ (t))^ - {X^ (t + r), X^ (t))^ 



, so that the squeezing spectrum 



^ can be written as 



ST (^) = 4 + 7out [SViy {u, - Vq v)] , 



(9) 



where the notation is self-explicative, see Eq. ((7j). This can be easily generalized to any 
pair of s-ordered two-time correlations. Taking into account that the s-ordered two-time 
correlation is nothing but 



(a (t + r),a(t)), 



^ (a (t + r) , a (t))^ + i-^ {a {t + r), a {t)), 



(10) 



with s G [—1,1], it follows that 
{a {t + r) ,a{t))j, - 



^ '-{a (t + r), a (t)), + {a (t + r), a (t)),, (11) 



s — s s — s 

for any s 7^ s'. Notice that Eq. (jHI) is retrieved from Eq. (jllj) when s = (symmetric 
ordering, which is obtained with the W distribution) and s' = — 1 (antinormal ordering, 
which is obtained with the Q distribution). Now, following the same arguments that lead 
to Eq. Q one gets 



1 — s' 1 — s 

-Vs {uj, V?) + Vs' {uj, ip) 



(12) 



s — s s' — s 

We see that the use of the V distribution is equivalent to the combined use of a pair of 
s-ordered distributions. 

The interest of this approach is that Eqs. (|7j) and (fT^ provide a way for comparing the 
predictions of a pair of s— ordered distributions, which we denote by Wg (Wi = P, Wq = W, 
W^i = Q), with that of the V distribution. This is interesting because the equation of 
evolution for a particular Ws needs not be of the Fokker-Planck t ype . For example, in 
the case we treat along this article (dispersive optical bistability [ll.ll^. llst). the equation 
of Ws includes additional terms (namely, third order derivatives) but for s = ±1, and. 



12, hj, but for s = 0. But 



in general, the diffusion matrix is not positive semidefinite 
tliese limitations do not necessarily prevent the use of these distributions as under some 
reasonable approximations their equations of evolution can be approximated to a Fokker- 
Planck equation (by neglecting the higher order derivatives in the Wigner case jl^ . llS^ or 
by limiting the study to a parameter domain where the diffusion matrix is well behaved in 
the Husimi case The point is that after making these approximations, Eq. ()12|) should 

provide not an exact but an approximate result. Then, by comparing the predictions of Eq. 
(fT^ to that of Eq. ((Tj) one could evaluate the influence of these approximations. 

In this article we shall make use of these approximations for the special case of dispersive 



optical bistability 
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12L llS^ . We then derive the fluctuation spectra from the linearized 
Langevin equations coming from s— ordered quasiprobability distributions, Wg, and compute 
the (linearized) squeezing spectrum. The main result we obtain is that, although any Ws 
obeys an approximate Fokker-Planck equation, and thus approximated Langevin equations 
can be obtained, the linearized squeezing spectrum given by Eq. (|T2|) is identical to that 
given by Eq. In other words, the approximations made in deriving Langevin equations 
from the approximated equations of evolution do not manifest in the linearized fluctuations 
spectra. We show further that the predictions for the (linearized) intracavity squeezing from 
any Ws is also exact. 

II. MODEL FOR DISPERSIVE OPTICAL BISTABILITY 
A. Master equation 

We shall adopt the model for dispersive optical bistability studied by Drummond and 
Walls consisting of a single-ended optical cavity containing a purely dispersive and 
isotropic x^^^ medium and pumped by a coherent field of frequency cu close to that of a 
cavity mode, uJc- The system Hamiltonian in the interaction picture reads 



H = h 



{6 — g) a) a + iEq {a) — a) — ^d^^a^ , (13) 



where Eq is proportional to the amplitude of the injected field, 9 = Uc — uj is a. detuning, and 
g = SeqItioIx/ (^^^) is the coupling constant, with V the quantization volume, e the medium 
dielectric constant and x = Xldi (^ = 1^2, 

we used a symmetrized Hamiltonian and this is the reason why the detuning is not 9 but 
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{6 — g): Had it been calculated in, say, normal or antinormal order, the detuning would have 
been 6 and {6 — 2g), respectively. Then, the correction g to the detuning is nothing but the 
modification of the cavity frequency ujc due to vacuum fluctuations as described with the 
different ordering choices. 

The intracavity fleld mode exits the cavity through the output mirror. Assuming weak 
coupling between the fleld mode and the rest of vacuum modes, which are treated as a 
reservoir (see, e.g., the master equation of the system at zero temperature reads 

p = Eq (a^p — pa^ + pa — ap) — i{9 — g) {a)ap — pa) a) (14) 
+ i— (a^'^a^p — pa^'^a^^ + — (2apa^ — pa^a — a^ap^ , 

where 7 represents the cavity losses of the fleld intensity and, as the cavity is single-ended, 
7out = 7- 



B. Quasiprobability distributions 

The equation of evolution for the s-ordered quasiprobability distribution for dispersive 
optical bistability was flrst derived by Vogel and Risken With our notation 



d_ 
dt 



Ws {a, a* 



92 



2 ^ daidaj 



D 



where ai = a, a2 = a* , and 



(15) 





= -£^0 




= Eq 







'^ + i{6 — sg) a + iga^a*, 
— i{6 — sg) a* — iga (a*)' 



9 1 — q 

tsga ~7 

-isg{a*f 



= D^f.) Making use of Eq. (HH) we obtain 



l-s' 



s — s' 

We shall make use of this property later. 



+ 



1 - s 
s' — s 



(16) 
(17) 

(18) 



(19) 
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For s = 1,0,-1, Ws corresponds to the P (Glauber-Sudarshan), W (Wigner), and Q 
(Husimi) distributions respectively as commented. An alternative quasiprobability distribu- 
tion to Ws is the so-called generalized P distribution , which we have already denoted by 
V. Its equation of evolution has been derived by Drummond and Walls [11]. It is given by 
Eq. (jl5j) with s = 1 after changing a* by the complex variable /?, which is independent of 
a and verifies {a) = {(3)* S. Hence in this representation the phase space is doubled with 
respect to the Ws representation. 



C. Fokker-Planck equation for the Wg distribution 



By construction all three V, P, and Q quasiprobability distributions formally obey a 
Fokker-Planck equation as they do not contain derivatives of order higher than 2. 

This is not the case for any Ws with s ^ ±1. We show next that, in spite of this fact, any Ws 
distribution verifies an approximate Fokker-Planck equation. This statement is equivalent 
to saying that, in some limit, the third order derivatives in E a. (11511 . existing unless s = ±1, 



can be neglected. For that we make a system size expansion 
very large value attained by the mean number of intracavity photons (a^a) 
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13|, which is based on the 
7/ 1^1 



. For example, by taking V = Icm:^, \x\ = 5-10 '^m Volt , e = Aeq, oj^ = 3-10 s 



one obtains \g\ ~ 10~^s~"'^ and, taking 7 = lO^s"""^, one has 7/ \g\ ~ 10^^. Now, normalizing 
time to 7 and a to \f^ijg-, one obtains an equation equivalent to Eq. (fTH|) in which the third 
order derivatives are multiplied by {gllf' ■, whilst the second order derivatives are multiplied 
hy g/'-y and the first order derivatives are of order one. Then, the neglection of the third order 
derivatives looks like a very accurate approximation. Notice, however, that the predictions 
of such a truncated equation may differ significantly from their correct values, as it occurs 
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13|. Finally note tha t, g iven the smallest value of g^ all A- 



with the tunneling times 
[Eqs. ()l(ipi7|l ] can be approximated by A, 
the system size expansion has been performed Eq. (jlSp becomes 



(s=0) 



; tna t, g i 

HQ, 



which we denote just by Ai. Once 



is) 



2^daidaj 

2„* 



Ws 



Ai = Eq — + iO^ a + iga^a 

— iO^ a* — iga («*)" 



Ao = En 



(20) 

(21) 
(22) 
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is a Fokker-Planck equation. However it is in fact a pseudo 
12 . llS^ as the diffusion matrix D*^*^ is not positive semidefinite in 



ai = a, a2 = a*, and D*^'') is given by Eq. (fTH|) . The symbol ~ is used instead of the equahty 
symbol in order to stress that Eq. (f^n|) is approximate (but for s = ±1). 

As announced Eq. ( 
Fokker-Plank equation 
general and then the equation cannot be interpreted as describing a generalized Brownian 
motion. That is not positive semidefinite is easy to see by writing Eq. (j20|l in terms of 
the real variables x = Kea, y = Ima. For these new variables an equation similar to Eq. 
pnjl is obtained with a diffusion matrix ©iy given by 



xy 



whose eigenvalues Sf^ read 



4" ~7±^ + 



The positive semidefiniteness of D^*) then requires that d^l^ > 0: 



< (23) 

(Remind that \af = x'^ + y'^.) Clearly, only for s = (Wigner distribution) condition 
fl23|) is fulfilled for any a. On the other hand for s = +1 (P distribution) condition ()23|) 
is never satisfied. In general, for any s ^ 0,1 that condition is verified inside a bounded 
region of the phase space and thus, D*^*-* is never, strictly speaking, positive semidefinite. 
Nevertheless if a is replaced by its classical steady value a, what is done for calculating 
linearized spectra as we do here (see below), D*^*) will be positive semidefinite whenever 
condition (jSHI) holds when applied to the classical steady state. This restricted condition is 
in fact verified in a bounded region of the parameter space. We note that this approximation 
was done for the case of second-harmonic generation by Savage It can be understood 
in the sense that one assumes that Ws is peaked around the steady state value and that 
the parameters of the system are such that a negligible part of the distribution violates the 
condition d^^^ > 0. Under this approximation 3^^^ is a well behaved diffusion matrix and 
the equation of evolution of Wg is a true Fokker-Planck equation for any s (7^ +1). 

The requirement of positive semidefiniteness of the diffusion matrix D'^'^^ comes from the 
fact that a and a* are complex-conjugate variables, and the noise terms in the Langevin 
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equations equivalent to the Fokker-Planck equation (see below) will not be complex- 
conjugated if this requirement is not fulfilled. Notice that this requirement is lifted in 
the case of the V distribution as a and (3 are not complex-conjugate variables, but in their 
mean. 



III. LINEARIZED FLUCTUATIONS SPECTRA AND SQUEEZING 

The Ito stochastic differential equations that are equivalent to the Fokker-Planck equation 
()20|) . assuming the positive semidefiniteness of the diffusion matrix, are 

^ ^ A (a) + («) • ^ (t) , (24) 

where a = (a, a*)^ , A = (Ai, ^2)^, ©("^ (a) = B^'^ (a) [B^") (a)]^, and the white Gaussian 
noise term £ = (^1,^2)"^ verifies {iiif)) = 0, and (^^ (t) (t')> = 5ij5{t-t'). (We note 
that cx should contain the label s, e.g. ct^'^\ as that stochastic variable is representation 
dependent. We avoid this labeling in order to not overburden the notation.) 

In this article we shall limit ourselves to the study of fiuctuations around the classical 
steady state a = a (the solution to A (a) = 0) in the linear approximation. The linearized 
Langevin equations read 

j^6a ~ A ■ (5q: + B(^) ■ ^ (t) , (25) 
, B(^) ^ B(^) (« = a) , (26) 

ai = a and 02 = a*. 

We are concerned with the calculation of the spectral matrix of fiuctuations, S*^*-* (to), of 
elements 

/+00 
dTe-'--{adt + T),a,{t))^, (27) 

(not to be confused with the squeezing spectrum) which, from the linearized Langevin equa- 
tions can be easily obtained by making use of Q| 



approx 



(uj) , (28) 
§i;)p_ (u) ^ (A + zul) ' (A^ - zul) , (29) 

where I denotes the 2x2 identity matrix and ID*^*^ = ©"^^^ (a = a.) and again the symbol 
~ stresses that the result is an approximation. The result for the different s— orderings is 
given in Appendix A. 



We remind that the Fokker-Planck equation obeyed by V is given by Eq. ()20j) for s = 1 
with the replacement a2 = a* (3. Then the full and linearized Langevin equations in the V 
representation are given by Eqs. (j24p and (|25p respectively, under the previous replacement, 
and the spectral matrix corresponding to the V distribution, S'^, is given by Eq. (j29|) with 
s = 1 (note that ]3 = a*): 

§^ {u) = (A + tul) (A^ - tut) . (30) 

We also remind that Eq. (j2Up for s = 1 is exact as the original Eq. (|15|) does not contain 
third order derivatives in this case; hence all symbols ~ must be replaced by = in this case, 
as in Eq. (jHUj) . On the other hand the normally ordered spectral matrix of fluctuations 
deflned as 




{d{t + T) ,a{t)) (a(t + r),at(t)) 
id^t + r) ,d{t)) (at(t + r),at(t)) 

equals, by definition, (cj). Thus Eq. (j!?n|) yields the exact normally ordered spectral 
matrix of fiuctuations, 

: § (cj) : = §^ (u) . (31) 
We now recall property (|TT|l that, together with definition (j2Z|), allows to state that 

§^ (,;) = (,;) + 1^§(^') (ou) , (32) 

s — s' s — s 

which, making use of Eq. (j2HI), can be approximated as 

^'^ (^) ^ ^^itox (^) + 737^— (^) ■ (33) 

Now, substituting Eq. (jSHI) into (jHHjl . and recalling property (fTTI|l and Eq. (jHOI), we observe 
that the approximate equality (jHHj) is a true equahty indeed. 

This means that, although the used Langevin equations come from approximated (trun- 
cated in general) Fokker-Planck equations, the spectral matrices Sipprox {00) given by Eq. 
(|29p provide the correct result. In other words, the third order derivatives present in the 
original pseudo Fokker-Planck equation seem to play no role on the correlations be- 
tween fiuctuations in a linearized theory. Moreover, the approximated Langevin equations 
()24|) ignore that the diffusion matrices can be non positive semidefinite, as discussed 
above. Nevertheless, as relation is not approximate but exact, and it holds for any 
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parameter set, even where D^^^ is not positive semidefinite, we conclude that the positive 
semidefinite condition on ID)^^^ is irrelevant in the following sense. We recall that the same 
occurs with the Glauber-Sudarshan P distribution: even if the pseudo Fokker-Planck equa- 
tion governing its evolution has a non positive semidefinite diffusion matrix one can 
nevertheless write down a corresponding Langevin equation, which yields the correct re- 
sult for the spectral matrix S [8]. The explanation for this was given by Drummond by 



introducing the generalized P representation j8[, which operationally amounts to substitute 
the complex conjugated variable a* in the pseudo Fokker-Planck equation verified by the 
Glauber-Sudarshan P by an independent complex variable (3. This suggests that our result 
can be understood in terms of "generalized Wg distributions", call them Wg- Should we 
substitute a* by an independent complex variable (3 in the original pseudo Fokker-Planck 
equation (fTK|l a positive semidefinite diffusion matrix would be not needed in order to derive 
corresponding Langevin equations (once the third order derivatives had been neglected). 
These Langevin equations for Ws would read as those for Ws but with a* — > (3, and the final 
expression for the spectral matrix, which would be exact in this generalized representation, 
would be given by Eq. (j^^ in the linear approximation, just as it happens in our case. 
The possibility of defining "generalized Wg distributions" should be studied, probably by 
defining Wg in terms of V, as Wg is defined in terms of P 3 • leave this discussion open 
as it is out of the scope of the present work. 

The above discussion implies that the squeezing spectra of the output field, given by Eqs. 
^ and ((121) with 

Vv (a;, ^) = \ [S^^ (co) e''^^ + (co) e^''^ + S^, (co) + S^, (co)] , (34) 
V. {u, ^) = i [SS (u) e-'^^ + Sii (u) + 5} J (u) + S^f {u)] , (35) 

are, obviously, the same. We shall not analyze here the properties of this squeezing spectrum 
as this analysis can be found in P]. We just quote in Appendix B the expression for Vg {on, (p). 
Nevertheless we want to make a comment on the amount of squeezing attainable inside the 
nonlinear cavity. This can be calculated by integrating the spectrum of fluctuations of the 
field quadratures, i.e. 

Vs{f) = 7r / dujVg{uj,<f). (36) 
We must take into account how the different Vg are related with the squeezing V {(f) = 
X^(t),X^(t)^ where is defined as X°^*', Eq. but replacing (^aout, with 
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(a, a^) . Performing the calculation one easily obtains, with the help of the commutator 
[a (t) , (t)] = 1 the result 

Vicp) = Vsicp)+'-. (37) 

The point is that, as any Vs Eq- (j35|) . can be computed from Sipprox (uj), Eq. p8|l . 

and Sipprox (uj) yields the correct result, any single Wg is useful for computing the intracavity 
squeezing. We note that V (ip) is minimum for a particular ip and is Vmin = |- This is a well 
known result: The maximum degree of squeezing attainable inside a nonlinear cavity, which 
happens at the bifurcation points, is half that of a coherent state. 

IV. CONCLUSIONS 

In this article we have discussed how the spectrum of squeezing of the field outgoing a 
nonlinear cavity can be derived from a combination of the spectra of intracavity fluctuations 
obtained from Langevin equations derived from Ws distributions. We have illustrated this 
for the special case of dispersive optical bistability. The interesting result is that the lin- 
earized spectrum of squeezing obtained is this way is exact in spite of the fact that no Ws 
quasiprobability distribution verifies Fokker-Planck equations but only approximate ones. 
We have also shown that the predictions for the (linearized) squeezing attainable inside 
the nonlinear cavity is correct when calculated with any Wg. The conclusion is that the 
linearized Langevin equations corresponding to a Ws representation are correct or, in other 
words, that the approximations made for converting the equation of evolution of Ws into 
Fokker-Planck equations do not manifest in the linearized theory, even if the diffusion ma- 
trix of the Fokker-Planck equation is not positive semidefinite. The latter has allowed us 
to conjecture the definition of "generalized Wg distributions", following the spirit of the 
generalized P distributions of Drummond. Of course, when going to the nonlinear regime, 
as in the calculation of e.g. tunneling times, one must be cautious about the truncation of 
pseudo Fokker-Planck equations containing third (or higher order) derivatives. 
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VI. APPENDIX A 



In this Appendix the expression for the spectral matrix Slpprox (oo) defined in Eq. (PU)) 
is given for any s— ordered quasidistribution Wg. The calculation needs the computation of 
the classical steady state (a = a,a* = a*), which is given by Ai = A2 = 0, Eqs. (j7H) and 
(jH}. After introducing the quantities 

3 - 
A = ^, /i = \g\ El ^le^f = ' a, = sign (g) , (38) 

the classical steady state is given by 

// = /[! + (/- Af] , (39) 

- ^ + (40) 
'l + il-Af 

The characteristic / vs. /i displays bistable behaviour for A > ^/S as is well known, and the 
values of the intensity / at the turning points of the characteristic, I = I±, are given by 

2A ± VA2 - 3 



(41) 



For I_ < I < the steady state is unstable; otherwise it is linearly stable. As the state 
equation implies that the pump power /i is univocally determined by / one can use the 
latter as the control parameter, and this is more convenient mathematically. 
Making use of Eq. one readily obtains {Q = 2uj /'-/): 

,^ 2 2iA-2IH3^^f<Sl^, (42) 



sS (^) 



{uj) = S[i i-u) , (44) 



where 



J =3(1 -/+)(/-/_). (45) 

Note that at the turning points of the characteristic (/ = /+ or / = /_) X = and all 
S^f (uj) diverge at = 0. 
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VII. APPENDIX B 



The expression for Vs {uj,ip), Eq. making use of Eqs. in Appendix A, 

reads 



where if) = 2 {(j) — ip) . 

As discussed (cu, ip) for s = +1 coincides with the corresponding expression calculated 
in the V representation: 

^, , , 2/(A-2/)cos^-r7/[2 + (l]2_j)]sinV' + 2/2 . . 

7Vp iuj, if) = — ^ n . (47) 

Note that this quantity is just 5'°"* {u) — \, Eq. ((Tj) (remind that 7o^t = 7)- 
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